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Abstract 

We evaluate the chiral condensate and Polyakov loop in two-dimensional QED 
with a fermion of an arbitrary mass (m). We find discontinuous m dependence in 
the chiral condensate and anomalous temperature dependence in Polyakov loops 
when the vacuum angle #~7r and m=0(e). These nonperturbative phenomena are 
due to the bifurcation process in the solutions to the vacuum eigenvalue equation. 
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The Schwinger model, QED in two dimensions has been a preferred theoretical lab- 
oratory for the study of physical phenomena such as chiral symmetry, gauge symmetry, 
anomalies, and confinement. JTJ-[|7J In a nontrivial topology it allows us to inquire about 
finite volume and temperature effects while keeping the computations infrared safe. ||- 



24[ Results at finite temperature (T) can be obtained by Wick rotating the solution on 
a circle S 1 of circumference L and replacing L by T . 

The theory is exactly solvable with massless fermions, but not with massive fermions. 
The effect of a small fermion mass (m/e <C 1) in the one flavor case is minor other than 
necessiating the 9 vacuum. [fl|, |6|, |I0[| The opposite limit of weak coupling, or heavy fermions, 
has been analized by Coleman. [J In this work we investigate physical quantities such 
as chiral condensate and Polyakov loop with no restriction on values of the parameters 
of the system. The effect of nonvanishing fermion masses has been investigated in lattice 
gauge theory and light cone quantization methods as well. P^HPTf 

The Lagrangian of the system is given by 

C=-\ V" + ^ 7 M (^ - eA^ - m(M + Aft) 

M = ^i(l- 7 5 )^, (1) 

where 7 M = (a±, io<i) and ip^ = We study the model on a circle of circumference 

L and boundary conditions 

A^t,x + L) = A„(t,x) 

iP a (t 7 x + L) = -ip a (t,x) . (2) 
The only physical degree of freedom associated with gauge fields is the Wilson line phase 

e w (t): 1,0,0 

e ie w (t) =expjie^ dxA^tix)} . (3) 

In the Matsubara formalism of finite temperature field theory boson and fermion fields 
are periodic and anti-periodic in imaginary time (r), respectively. Mathematically, the 
model at finite temperature T = is obtained from the model defined on a circle by 
Wick rotation and replacement L — >• /3, it — ► x and x «-> r. The Polyakov loop of a charge 



q in the finite temperature theory corresponds to the Wilson line phase: 

P g (x) = exp fa J Q drA (T, x)\ exp ji^ W (*)| • (4) 
We bosonize the fermion in the Coulomb gauge in the interaction picture defined by 



a massless fermion: [110, |T8 



Jj U % \ _ \ Q e ±i{q±+2-Kp±(t±x)/L} . e ±n/I^±(t,z) . 
V 7 ! 



C + = 1, , C_ = exp{?7r(> + - p_)} , 

oo 

= E(4™) _1/2 {c ± ,„e- 2 ™^/ L + h.c.} , 

n=l 

e 2« P± |ph ys ) = |phys) , (5) 

where [<7±,p±] = i, and [c± in , c4,m] = ^nm- The : : in (|]) indicates normal ordering with 
respect to (c n ,cjj. In physical states p± takes an integer eigenvalue. 
The Hamiltonian in the Schrodinger picture becomes 

7T 

Htot — H, o + + H mass — — 

bL 



H^= j dx - : 



£ 1 r P 2 

n 2 + (0') 2 + -0 2 



7T 



^mass = f L dxm(M + M^) . (6) 







The conjugate pairs are {p, g} = {±(p+ + p_), q + + {p, q} = {p + -p_, ±(g + - q_)}, 
{Pw,@w}) an d {n, = + + Note that (0,11) fields are subject to conditions 

Jq dx(p(x) = = f L dx U(x). The mass operator is given by 

M = -ClC+ ■ e ~ 2niiix/L e iq ■ L _1 iVo[e iV5F<s ] (7) 

where iV^[- • •] indicates that the operator inside [ ] is normal-ordered with respect to a 
mass pi. 

As [p, -fftot] = 0, we may restrict ourselves to states with p = 0. p takes integer 
eigenvalues in this subspace. The Hamiltonian (Bp posseses a residual gauge symmetry 



0w - ► @w + 27r and p — > p — 1 generated by Z7 defined aspHj] 



17 = exp (2tuP w + ig) , [17, H tot ] = . (8) 

The ground state is the 9 vacumm: @-@ 

f/|$vac(^)) = e i9 |$ vac (^)) . (9) 

To determine the vacumm wave function we must solve the eigenvalue equation 

(H + # mass )|$ vac (#) ) = E\$ vac (6) ) . (10) 

The fermion mass term i7 mass changes the mass of the boson field <fi from \x = g/^/ti to 
/ij. The vacuum is defined with respect to the physical boson mass H\. Making use of 



10, 13 



No [e 47r< ^] = B(jiL) N^e 1 ^] 

B{z) = ± exp L + 1 - 2 f 1 , (11) 

the mass operator in fl7|) is accordingly written as 

M = -ClC + ■ e - 27Tipx/L e ig ■ L- l B(^L)N^ [j^] . (12) 
With this understanding we write the vacuum wave function, taking (||) into account, 

as 

|$ vac (£)> = -±=£ / dp w \p W7 n)e- ine+2 ™™ f(p w ) 
v 2tt n J 

J dp w \f{p w )\ 2 = l . (13) 

Here \pw, n ) is an eigenstate of P\y and p. Since ( p' w , n'\e ±iq \pw, n ) = 8(p' w — pw) 8 n ',n±i, 
the vacuum eigenvalue equation ( |T0| ) is reduced to a Schrodinger equation 

d 2 



dp 



w 



, +V(pw)\ f( Pw ) = ef(p w ) (14) 



where 



V(Pw) = u 2 Pw ~ K cos(6 — 2np w ) 
4 



e 2 

uj = nfiL , k = 4:7imLB({iiL) , /i 2 = — (15) 

7T 



and e = 27iLE va _ c + ^n 2 . 

To determine the boson mass /ii, we expand H mass in (|6|) in power series in 0. In the 
vacuum 

" L , 2nm 



H mass - / dx--B^ 1 L)(e^ + e-^\ 
Jo L 



i.2 



rfx r VP ; ( cos (9 - 2np w )) f 2 16 

lo L 1 

where (F(p w )) f = J dp w F(p w ) \f(p w )\ 2 - Hence 

/i* = /i 2 + J^L^ ( cos(0 - 2 W ) ) y . (17) 

(|l^) and (|17D must be solved simultaneously. We have a Schrodinger problem in which the 
potential needs to be determined selfconsistently. Wave functions f(pw) for typical 
values for T//x, rn//i, and 9 are displayed in fig. 1. 
The chiral condensate is given byQ 

= -2L- l B{^L) ( cos(9 - 2np w ) ) f . (18) 
Combining (O) and (ITH), one finds 

fij-fi 2 = -4vrm(^) e , (19) 

which is a PCAC relation. 

The Polyakov loop is, from (^), 

/ p \ -I J(q/e)®w \ 



q 

for - ^ an integer 



dpw f(Pw)*f(Pw ~ -) for - = an integer. 

-oo e e 



(20) 



3 We have defined the mass operator M by (Q), independent of a mass m. Consequently (ipip)^ 
even in the limit e^O. It may be appropriate to define a physical M by M phys =Af — ( M ) m e=0 so that 
(^^)^| lys vanishes in the free theory. The composite operator M phys thus defined depends on m. 
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Figure 1: Wave functions \f{pw)\ 2 - ( a ) vnj n— 10, 1, and 0.1 corresponds to k — 3.61 • 
10 _1 , 3.19 ■ 10~ 2 , and 3.18 ■ 10 -3 . (b) 9 = n~ indicates that the 9 value is less than, but 
is very close to ir. At m//x=l, 9=0 and 7r corresponds to k— 2.24 and 2.19. At m//x=10, 
k=165 for both #=0 and 7T. 



The potential, V(pw), in Eq. (0) consists of two terms; harmonic oscillator and cosine 
potentials. The strength of the cosine term, k, is 

_ fe 7 m/iiL 2 for [i\L > 1 
K ~ XatiuiL for niL < 1. 

Depending on relative strength, the behavior of the ground state wave function is quite 
different. If uj 2 ^> k, the potential is approximated by the harmonic term, ie. f{pw) — 
(u/'K) 1 / 4 e~ ujp vi-/ 2 . The condition is satisfied if raj\i <C 1 <C \iL or if m/fi <C fiL <C 1. In 
this regime ( cos(9 — 2irp w ) ) j = e~ n ^ L cos 9 so that 

[i\ = \J + (me 7 cos 9) 2 + me 1 cos 9 

6 



— e 7 
(i/ji/j) g = ——//i cos6» 

(Pe> fl = e-^ /4T ~ for - « 1 « ftL = % 

a T 



(22) 



and 



£ = + —^ e -*/*z> cos 9 

= --e-^ cos 9 
(P e )e = e-^ /4T ~l for-« / ,L = ^«l , (23) 

In the opposite limit k ^> uj 2 , the cosine term dominates. However, the harmonic 
potential cannot be ignored as it lifts the degeneracy of the cosine potential. In the large 
volume limit the harmonic potential selects one of the minima of the cosine potential at 
p w = 9/2-k where 9 = 0- 2n[(9 + n)/2n]. Hence V ~ \e l m^ l L 2 (2-np w - 9) 2 so that 
f{p w ) = ^ L y/^ e -i^L(2np w -Sr/8^ where ^2 = 2m/Ul e 7 . [For 9 ~ 7r (mod 2vr), / has two 

peaks at pw ~ This leads to ( cos(# — 2npw) )/ = e~ n '^ L ~ 1. Consequently 



Hi = jl = 2e 1 m 

e 27 

(V'V'/fl = m 



TT 



( P e >e = e^ e7m / 2T for m » /i , //L = | > 1 . (24) 

Notice that the chiral condensate increases linearly with m. However, there is no 9 
dependence to the leading order. (See the previous footnote, too.) 

At high temperatures (T = L~ x ^> /i), lo 2 j k = 7r/z 2 /4mT. So long as m ^ 0, eventually 
the cosine term dominates. However, the both terms in the potential become small in 
this limit. A good estimate is obtained by treating the cosine term as a perturbation. 
Numerical evaluation supports the result. (See fig. 3 below.) 

Write Eq. (|T4]) as (H + V')f = ef where V is the cosine term. Eigenstates of H are 
denoted by {|n )} (n = 0, 1, 2, • ■ ■). Then 



( n| cos (9 - 2np w )\0 ) = (<** + ( l)"^*) ^ 1 , 2* n * 

^ vn! 
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Figure 2: m//i dependence of chiral condensates at given T//i and 9. 



f is given by \f ) = |0 ) — J2 n =i \ n ) ( ti| V|0 )/2un so that 
27rp w ))/ = e^ 2/<Jj cosfl 



cos 



-2tt 2 /u 



2n 2 /u> e Z -1 



dz h cos 26 I dz ■ 

z Jo z 



(26) 



2uj [ Jo 

For T = L^ 1 <C \l (ui ^> 1) the first term in (^6|) dominates over the rest to reproduce 
@. For T > // (a; < 1), Jo 2 "^ ^ (e z -l)/^ ~ (^/2vr 2 )e 2 " 2 / w so that ( cos(6-2tt Pw ) )/ ~ 
m/irT. We obtain 



fil = /j 2 + 8m 2 

<**>, = -- 

7T 



(P e ) e = e-^ 4T 1- 



6.60 m 



-7rT/ M 



cos^ I for T ^> /i, m. 



(27) 



In the intermediate range of parameter values Eqs. ( [14]) and ( p|) must be solved 
numerically. The computational algorithm is the following. With given (fiL,m/ fi,0), we 
first assign an input value for k = n in . The potential in (|14|) is specified with (fiL, n in , 6). 
Eq. ( [14] ) determines f(pw), from which one can determine /xi, solving (0). With this 
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Figure 3: T dependence of chiral condensates at fixed m/fi and 9. 



new Hi one recomputes k = n out by (|15D . Schematically 



V(pw) -> /(pw) 



^out • 



(28) 



/t out must coinsides This gives a consistency condition to determine k with given 
(p,L,m/fi,9). 

In fig. 2 chiral condensates are plotted as functions of m/fi. Asymptotically (m ^> fi) 



the behavior is given by (24|). T dependence of chiral condensates is displayed in fig. 3. 

. The T dependence of 



24| Low and high temperature limits agree with 



and 



chiral condensates with given m and 9 is smooth. This is consistent with the Mermin- 
Wagner theorem that in a one-dimensional system there is no phase transition at finite 
temperature. |28 



T dependence of Polyakov loops is displayed with various values of mj \x in fig. 4. Notice 
that at 9 = the curve smoothly changes as m/fi. However, around 9 = ir nontrivial T 
dependence is observed for m > //. The origin of this behavior is traced back to the wave 
function f(pw) m the corresponding problem on S 1 with L = T _1 . f(pw) at moderately 
low temperature has two dominant peaks at pw = ±(1 — e) around 9 = it, whereas it 
has only one dominant peak at pw = for 9 = 0. As (JZHp shows, the Polyakov loop is 
determined by the overlap of f(pw)* an d fipw ~ !)• Hence it vanishes quickly for 6 1 = 0, 
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Figure 4: T//i dependence of Polyakov loops with given mj\x and 9. Two-step behavior 
is observed at 9 ~ 7T. 



but approaches ~ .5 for ~ 7r at low, but moderate T . As T further gets lowered, the 
two peaks become narrower and sharper. When the width of the peaks becomes smaller 
than e, the overlap of the wave functions and Polyakov loop vanish. If 9 is not exactly 
tt, but is very close to tt, the asymmetry in the potential, enhanced by the factor m/T 2 , 
becomes important at sufficiently low T and the wave function has a sharp peak around 
one of the minima. See the wave functions displayed in fig. 1. In the numerical evaluation 
presented in fig. 4, the computer picks a value for 9 which is not exactly tt. The transition 
from the plateau (~ .5) to zero at 9 ~ 7r is caused by this change of the wave function. 
This explains the two-step behavior observed in fig. 4. 

In fig. 3 we observe that the pattern of the T dependence of (ipip) changes as m/fi 
is increased. At 9 = tt, (ipip)/fi decreases (increases) as T//i increases for m//i = .1 
[mjii = 1). 

Indeed there arises a discontinuity in the m dependence of chiral condensates at low 
temperature at 9 = tt, at least in our approximation scheme. We have displayed it at 
T j ii = 0.03 in fig. 5. We observe that (ipip)/^ discontinuosly changes at m/fi = .437. 

To understand the origin of the discontinuity, we consider a solution = K out in 
(p8|). Write K out = g(K in ). We are looking for a solution to g^n) = k, or a fixed point 
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Figure 5: A discontinuity in chiral condensates is observed at the m/jj, ~ .437 when 
T/n = 0.03 and 9 = ir. 

Kfi x of q(k). As m/ji varies with given T//i and 9, q(k), and therefore kr x change. At 
9 = it, there is a critical value for m/fj, at which the fixed point bifurcates. In a certain 
range of m/fi, there appear three fixed points, two stable and one unstable. Among the 
two stable fixed points, one of them has a larger chiral condensate and therefore a lower 
energy density, corresponding to the vacuum. This bifurcation induces a discontinuous 
change in chiral condensates as m/fi varies. We have displayed the mapping K out = g(K m ) 
in the critical region in fig. 6 (a). We observe that saddle node bifurcation takes place at 
m/n = 0.4368. 

At very low T the critical mass m c j ' ji can be determined analytically. Suppose that 
m = O(n). For /xL = /i/T > 1, k = e'm^L 2 . At 9 = n, 

V = (irfiLp w ) 2 + ncos2irp w . (29) 

There is always a solution which satisfies (fiL) 2 > 2k or ji 2 > 2e l m\X\. In this case f{pw) 
is sharply localized around pw = 0. This yields 



Hi = y fi 2 + (me 1 ) 2 — me 1 . (30) 

If (^L) 2 < 2k or fi 2 < 2e 7 m/ii, f(pw) is localized around ±pw 7^ where 

_ 2mii 1 e l . 

2np w = — sm2irp w 

H 2 
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Kin Kin 

Figure 6: Ki n -K on t plots, (a) At T/fj, = 0.03 and 9 — ir. As raj ji changes, the fixed point 
bifurcates at m/fi = 0.4368. (b) With a fixed m/fi = 0.4365, T/fi is varied. The curves 
remain almost critical. 



/ii = y /i 2 + (me 7 cos 27ipw) 2 — rne 1 cos 2tcpw (31) 

A solution to (|31|) exists only for m/fi > 0.435. This solution corresponds to a bigger 
fii or —(ipifj), and therefore to a lower energy density. In other words m c /fi = 0.435 at 
T = 0. Numerically we have found m c //i = 0.435,0.437, and 0.454 at T/fi = 0.01,0.03 
and 0.07, respectively. Above T/fi = 0.12 the function g(x) has only one fixed point for 
all values of m/fi so that the discontinuity disappears. 

It is not clear, however, if the discontinuity discovered above is real in the full theory, 
or just an artifact of the approximation in use. In determing the boson mass fi\ in ( |l~6|) 
and ([17]) we have ignored nonlinear terms in the field in H mass , retaining only the 2 
term. Those nonlinear terms are expected to affect the boson mass. 

Furthermore, the Mermin- Wagner theorem ensures that there should be no discon- 
tinuity in when T j [i varies with m/fi kept fixed. The discontinuity in the m/fi 
dependence is consistent with the Mermin- Wagner theorem only if m c /fi is universal, 
being independent of T//i. Numerically we have found that m c //i is almost universal, 
although there appears tiny dependence on T/fi. In fig. 6 (b) we have plotted n in -K out 
with a fixed m/fi = .4365 at various T/fi. The curve is critical at T/fi ~ 0.03, while it 
is slightly off at higher or lower temperature despite it may be very hard to see visually. 
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It is striking that the critical value m c / /i is almost insensitive to T//x to such a degree of 
accuracy. 

There are two possible senarios when the whole interactions are taken into account. 
It may turn out that the discontinuity in mj\x disappears, being replaced by a crossover 
transition. Or the discontinuity is real, taking place at a universal mass value m c / '/i. 
Further investigation is necessary to determine which picture is right. 

In this paper we have evaluated T- and m-dependence of chiral condensates and 
Polyakov loops. We have demonstrated that at #~7r there appears anomalous behav- 
ior when m/fj, is 0.4 ~ 0.5. Mathematically these anomalous phenomena are related to 
the bifurcation process in the solutions of the vacuum equation. Electromagnetic interac- 
tions and fermion mass collaborate to induce anomalous behavior of the vacuum solutions. 
This is reminessent of chaotic dynamics in nonlinear systems. More detailed analysis will 
be reported separately. 
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